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Abstract 

We discuss from a bi-Hamiltonian point of view the Hamilton-Jacobi separability of a few 
dynamical systems. They are shown to admit, in their natural phase space, a quasi-bi- 
Hamiltonian formulation of Pfaffian type. This property allows us to straightforwardly 
recover a set of separation variables for the corresponding Hamilton-Jacobi equation. 

1 Introduction 

The notion of quasi-bi-Hamiltonian (QBH) systems has been introduced quite recently Q; 
it originates from the study of dynamical systems which are associated with two compatible 
Poisson tensors (Po,Pi), but do not admit a bi-Hamiltonian formulation, in the given phase 
space, w.r.t. to these tensors. 

At least for a class of such systems (the so-called Pfaffian QBH) it has been shown Q how 
to construct a set of separation variables for the corresponding Hamilton-Jacobi equation, so 
that a Pfaffian QBH system can be integrated by quadratures. In our opinion, this relation 
between QBH formulation and separation of variables is remarkable. 

On the other hand, since at present there is not, at the best of our knowledge, a satisfactory 
general scheme encompassing these properties, we believe that it may be preliminarily useful 
to study some concrete examples of such systems; as a matter of fact, it turns out that the 
QBH formulation is shared by a few classical separable systems (as for other cases previously 
considered in the literature, see ^ Q ||, g]). In the next section we briefly review the main 
definitions and results to be used in Sect. 3, where the following systems are presented: the 
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Kepler problem with a homogeneous force, the Euler problem with two fixed centers and an 
elastic force, the motion on the ellipsoid in an elastic potential and a solvable n-body problem 
introduced in |Q. 



2 Bi-Hamiltonian manifolds and quasi— bi— Hamiltonian 

systems 

A bi-Hamiltonian (BH) manifold || is a smooth manifold M endowed with two compatible 
Poisson tensors Pq, Pi : T*M i-s- TM {TM and T*M being the tangent and cotangent bundle 
of M, respectively). As it is known, if M is even dimensional (dim M = 2n) and at least a 
Poisson tensor, say Pq, is invertible, then N = P\ P^ 1 is a Nijenhuis tensor || (a hereditary 
operator in the terminology of In particular, if N is maximal, i.e., it has n functionally 

independent eigenvalues (Ai, . . . , A„), one can introduce a set of canonical coordinates (A; fi) 
(A := (Ai, . . . , A„); fi := . . . , (J, n )) referred to as Darboux-Nijenhuis coordinates such that 
Po, Pi and N take the matrix form 



Po = 



I ' 


Pi = 


A 


N = 


' A 


" 


-I 


-A 





A 



(1) 



where / denotes the n x n identity matrix and A := diag(X%, . . . , A„). The first n coordinates 
are just the eigenvalues of N whereas the remaining ones can be constructed by quadratures 
UM- The above form will be referred to as the canonical form of the previous BH structure. 



2.1 Bi— Hamiltonian systems 

A vector field X is said to be bi-Hamiltonian w.r.t. a pair of Poisson tensors (Po, Pi) if there 
exist two smooth functions ho and hi such that 



X = P dhi = Pi dho , (2) 

d denoting the exterior derivative. 

If N exists and it is maximal then a BH vector field is completely integrable, a set of 
independent involutive integrals being just the eigenvalues (Ai,...,A„) Q. Conversely, a 
strong condition has to be satisfied by a completely integrable system in order to admit a BH 
formulation in a neighborhood of a Liouville torus, at least if one searches for a second Poisson 
tensor compatible with Po Hereafter, this kind of structures will be referred to as 

standard bi-Hamiltonian structures. 

Nevertheless, the property of Liouville integrability can be related with geometric formula- 
tions which are actually different from the BH one w.r.t. a standard BH structure. This can 
be done in (at least) three distinct ways. 

i) Searching for a BH structure (Qo, Qi) not including the Poisson tensor Po Jl^]; hereafter 
such structures will be referred to as alternative BH structures. 
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ii) Admitting a degenerate BH formulation; this is the case, for instance, of the rigid body 
with a fixed point (T^, [Tt], [l^] and of the stationary flows of the KdV hierarchy Jl^, [2(J . 

hi) Searching for a QBH formulation of the vector field X w.r.t. a standard BH structure 

Si- 



2.2 Quasi— bi— Hamiltonian systems 

In connection with the third approach, we recall that the vector field X is said to be quasi-bi- 
Hamiltonian w.r.t. a pair of Poisson tensors {Pq,P\) [Q if there exist three smooth functions 
H, K, p such that 

X = P dH = -Pi dK . (3) 
P 

One could say that the given dynamical system is Hamiltonian also w.r.t. the Poisson tensor 
Pi, provided that a nontrivial change in time: dt — > dr — (1/p) dt is introduced. In particular, 
a QBH vector field X is said to be Pfaffian if N exists and p — fTCLi Ajj e -' ^ is just the 
product of the eigenvalues of N. The relevance of QBH vector fields is based upon the following 
two results. 



Any completely integrable system with two degrees of freedom admits a QBH formulation 
in a neighborhood of a Liouville torus |Q . 

Any Pfaffian QBH vector field with n degrees of freedom is separable (in the sense of 
Hamilton-Jacobi) in the Darboux-Nijenhuis coordinates Indeed, the general solution 
of Eq. (|J), written in these coordinates, is 

g = 2^ rr" .n.-iA ■ K = 2^ rr» x a (ft : = (4) 



^ IIi=i (Ai-Aj) ' ^n"=i(Ai-A,) 



i=i 



where each function/^ is an arbitrary smooth function, depending at most on one pair of 
variables (Ajj /ij). (Obviously enough, if n — 1 it is: rij=i(Ai — Aj) := 1, = 1, Pi = XPq, 

H = K = /(A, /it)). A remarkable feature of H and if is that they are separable as they 
verify the Levi-Civita condition |2l) , so that the corresponding Hamilton equations are 
integrable by quadratures. We stress the fact that, owing to the arbitrariness of fi, the 
functions H and K (^) provide a class of separable functions generally different from the 
known Stackel class, quadratic in the momenta (e.g., see [^2], p. 101]). 

The results of 0] have been completed in jij , where it has been shown that a QBH vector field 
X admits n integrals of motion in involution Fk (k = 1, . . . , n) 

*"fe»>*iiy=i(^-^) ' 1 J 
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Ci , . . . , c n being the coefficients of the minimal polynomial of ./V 

n n 

in particular, F\ = —H,F n = (—l) n K. Furthermore, each function turns out to be separa- 
ble. 

Remark. Let H, K, p be of the form (Eh, i.e., the general solutions of the QBH Eq. (^) in 
the Pfaffian case. If <p , tp : R — ► R arc C 1 functions with nonvanishing derivatives ip', ip' , let 
us define on the phase space the functions <J>, VP given by 

$(A; /i) := (A; /*)) . *(A; m) == V>(#(A; M)) • (7) 
Then $, $ are solutions of the QBH Eq. (^) with a function p given by 

Hence, for each QBH system of Pfaffian type there is a class of QBH systems of non Pfaffian 
type, with p of the form Moreover, it is straightforward to show that $ and '5 are separable, 
since they satisfy the Levi-Civita condition, provided that this condition be fulfilled by H and 
K (as in the Pfaffian case), whatever the particular form of the functions <p and tp. 

Viceversa, let us search for the general solution of Eq. (||) w.r.t. H, K and p = pf'(H)/g'(K) 
(p = n"=i Ai), with / and g arbitrarily chosen functions; it is easy to prove that H and K are 
given by 

H = f-\H), K = g- 1 (K) (9) 

with H and K of the form (^) . 

This kind of generalization of the Pfaffian case has been recently considered in |23| only for 
the case n — 2. As it is evident, the above results hold for QBH dynamical systems with an 
arbitrary number of degrees of freedom. ■ 



2.3 The "origin" of some QBH vector fields 

Let us describe a possible situation in which some interesting QBH vector fields arise. Let 
(M,T'o,'Pi) be a BH manifold and X the vector field of a given dynamical system on M, 
Hamiltonian w.r.t. Vq. If neither Vo nor V\ is invertible, a possible way to analyse the inte- 
grability of X is to eliminate the Casimir functions of one Poisson tensor, say Vo, by fixing 
their values; of course, both Vo and X can be restricted to a symplectic leaf 5*0 of Po, so that 
X is still a Hamiltonian vector field on So. However, if P\ cannot be restricted to So, the BH 
formulation (||) is lost on So, even if X is BH on M. As a matter of fact, in a few cases the 
following situation occurs: 
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i) there is a fibration tt : M — » M' = M/7T such that both Vq and Pi are projectable 
along 7r; since 7r turns out to be transversal to So and Si (a symplectic leaf of Pi), the 
quotient space M' and the symplectic manifold So are diffeomorphic and So itself is a BH 
manifold, with Po and Pi invertible, Po and P\ denoting the reduced tensors of Po and V\ 
respectively. If this is the case, and if the eigenvalues of the Nijenhuis tensor N = PiP -1 
are independent, then one can introduce the Darboux-Nijenhuis coordinates on So. 

ii) There is a function p such that the restricted field X admits the QBH formulation (||), 
with H and K given by the restriction to So of integrals of motion of X. The interest of 
this result is that if X is Pfaffian, then in the Darboux-Nijenhuis chart it is separable ||. 

We remark that the situation described in i) and ii) is, so to say, "experimental" , a sound 
theoretical foundation of these results being lacking. The peculiarity of the above reduction is 
given by the fact that two different geometric processes are used simultaneously: the restriction 
for the vector field and the projection for the BH structure. Due to this fact, one maintains 
the BH structure but loses the BH formulation for the vector field, recovering in some cases the 
QBH formulation. This happens, for instance, for: 



the integrable Hcnon-Hcilcs system and its multidimensional generalizations obtained by 
reduction from the stationary flows of the KdV hierarchy Q ; 

a class of permutationally symmetric potentials recovered from the restricted flows of the 
coupled KdV systems, the most representative member being the Gamier system B. 



Both classes of dynamical systems live on a BH manifold M of maximal rank (dimM = 2n+ 1) 
and their QBH formulation is obtained by the reduction to a symplectic 2n dimensional manifold 
according to the above scheme i) , ii). The study of examples with BH structures of non- 
maximal rank such as the stationary flows of the Boussinesq hierarchy will appear elsewhere 
p^j. In this regard, we recall that the geometry of (2n + l)-dimensional BH manifolds of 
maximal rank has been completely described by Gelfand and Zakharevich ]26|| , whereas there 
is not a similar analysis for non-maximal BH manifolds. 

For the sake of clarity, we recall that we are considering from the very beginning a well 
specified vector field X which is associated to the dynamical system under investigation. So, 
the situation described above does not contradict the following known results: if one has been 
given a BH structure with Po invertible, one can construct a BH vector field whose Hamiltonian 
functions are the traces (or the eigenvalues) of the Nijenhuis tensor N |l2[ ; moreover, a QBH 
vector field can always be constructed for any maximal Nijenhuis tensor J25|. However, it turns 
out that these BH and QBH fields do not in general coincide with the vector field X of the 
given dynamical system. 



3 Examples of dynamical systems with QBH formula- 
tion 

We present four examples of integrable systems which can be given a QBH formulation of 
Pfaffian type w.r.t. a standard BH structure. 
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The first three systems (see, e.g., M, p. 126-129]), are defined on the cotangent bundles of 
Riemannian manifolds; therefore they are of Stackel type and the Nijenhuis tensor of their QBH 
formulation can be constructed by lifting, in a suitable way, a conformal Killing tensor |^| from 
the configuration space to the corresponding cotangent bundle. (Such a construction will be 
the subject of a further publication). The dynamical system considered in the last example 
is naturally Pfaffian QBH, the physical coordinates being just Darboux-Nijcnhuis coordinates 
w.r.t. a standard BH structure written in the canonical form (|l]). 

3.1 The Kepler problem with a homogeneous force (Lagrange 1766) 

Let us consider the classical problem of a particle in the plane under the influence of the Kepler 
potential and of a homogeneous field force. The Hamiltonian function is 



H 



/ 2 2 \ ^ 



bq2 



(10) 



where {q%, qi'iPx^P'i) are the cartesian coordinates of the particle and the conjugate momenta, 
respectively; a and b are real constants. The vector field is X — Po dH, Po being the canonical 
Poisson tensor. There is a second independent integral of motion 



v Pi, s Qg2 . bql 

K = ~7C\Q2Pl - HP2) - — F5 , + —T 
2 2y/qJ + qi 4 



(11) 



which allows us to give X a QBH formulation. Indeed, we can write X = -Pi dK with p = — 2^ 
and 












qi 








<7i 


2g 2 





-qi 





-pi 


-qi 


-2(72 


Pi 






Pl = 2 



Since the minimal polynomial of the Nijenhuis tensor N = P\Pq 1 is 

m(A) = A 2 — q2\ 



(12) 



(13) 



one easily checks that X is a Pfaffian QBH vector field. The Darboux-Nijenhuis coordinates 
can be constructed following |0|; they are 



Al 2 



\(q2T 



Ml, 2 = P2 



— W2 ± 

qi 



(14) 



One can easily verify that (Ai, A2) are just parabolic coordinates in the plane with focus at the 
point q\ = 92 = and axis the (72-axis. In the above coordinates the two integrals of motion 
read 
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H 



Ai/i 2 - 2b\j 



X 2 ^ 2 — 26A 



K = X 



2(Ai-A 2 ) 
Aim? - 26A^ 



+ A 



2(A 2 -A!) 
A 2 jul - 26A 2 + a 



2(Ai-A 2 ) 2(A 2 -Ai) ' 

On account of the general result proved in M, one recovers that H and if are separable. 



(15) 
(16) 



3.2 The Euler problem with the two fixed centers and an elastic force 
(Euler 1760-Lagrange 1766) 



The Hamiltonian function for this problem can be written as follows 



H 



1/ 2 , 2 )A b 2 b ? 

2 lPl P2> V^ + fe + c) 2 Vqf + (i2-cr 



(17) 



where: (q\, q 2 ; pi, p 2 ) are as in the previous example; the two centers are at the points F\ = 
(0, — c) and F 2 = (0, c); 61, b 2 , k are real constants. The corresponding vector field is X ~ Po dH. 
As above, also this system does not admit a BH formulation w.r.t. a second Poisson tensor Pi, 
but it can be given a QBH formulation. Indeed, one can write Eq. (^|) with K given by the 
integral of motion [a\ S R, a 2 = a\ + c 2 ) 



a 2 2 ai 2 1 2 b 1 (a 2 + cq 2 ) b 2 (a 2 ~cq 2 ) k 2 2 

K = —p 1 + —v 2 -^(<liP2-q2Pi) + , === + ^== = +-(a 2 g 1 +a 1 g 2 ) , (18) 

222 V?i + (<?2 + c) 2 V9i + (<?2 - c) 2 2 



p = (aia 2 — a 2 q\ — ai9 2 ) and the following Poisson tensor P\ 



Pi = 






-{ai-q\) 

9l92 



fli - q\ 

-9l?2 

9i 92 

-{ai-ql) 91P2-92P1 



-91 92 
«2 - 9 2 
-qiP2 + q2V\ 




The minimal polynomial of the Nijenhuis tensor N = P\P n 1 is 



(19) 



(20) 



m(A) = A 2 + (g 2 + q\ - 01 - a 2 )A + aia 2 - a 2 g 2 - aig 2 , 

so one easily recognizes the Pfaffian property of X. The Darboux-Nijenhuis coordinates are 
given by 



Ai j2 
Mi, 2 



(ax + a 2 - q\ - ql) \J {q\ + q\ - a\ - a 2 ) 2 - A{a x a 2 - a 2 q\ - a^qfj 



Pi 92 



2 

P29l 



49i9 2 

(Pi92 - P29l) 
4(a 2 - ai)qiq 2 



T 

(Pi92 ~P29i)(9i + 9 2 ) 
4(a 2 - ai)gig 2 



(a 2 - ai)(a 2 - 01 + 2g 2 - 2 g 2 ) + (9? + 9I) 2 



(21) 



(22) 
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One can verify that (Ai, A2) are just the Jacobi elliptic coordinates in the plane (already known 
to Euler) with foci at the points F\, F 2 and axis the g 2 -axis. In the above coordinates the two 
integrals of motion read 



H 



-2(ai - Ai)(a 2 - Ai)/^ - (h + b 2 )y/a 2 - Ai - |A| 
Ai — A2 

+ -2(ai - A 2 )(a 2 - Agj/jj + (b 2 - h) sqrta 2 - A 2 - |A| 



A 2 — Ai 

— 2Co,i — Ai Von — Ai 

K = X 2 



2(oi - Ai)(a 2 - Ai)^ - (h + b 2 )V^A7 - |Ag 



Ai — A 2 

-2(ai - A 2 )(a 2 - A 2 )m| + (62 - h)Va 2 - A 2 - |A| 

+ Ai (24) 

A 2 — Ai 

Also in this case they are separable, since they have the general form (J4J) . 

3.3 The motion on the ellipsoid in an elastic potential (Jacobi 1843) 

Let us consider a harmonic oscillator on the n-dimensional ellipsoid 



i=0 

where Xi(i = 0, 1, • • • , n) are Cartesian coordinates in R n+1 and ao < a\ < a 2 < • • • < a n 
are positive real constants. If yi are the conjugate momenta, the Hamiltonian function is 
H = I SiLo(?/i + x ? )• Let us cons ider the canonical transformation (a;,, J/i) 1— > (A;, /Ltj) associated 
to the point transformation (x^ t— > (Xi), where Ai are the generalized elliptic coordinates in 
R n+1 [^7j defined as the (n + 1) roots of the equation 

n 2 

E^ L T = 1 - ( 26 ) 

The ellipsoid (|2^) is given by the submanifold Ao = 0, and the Hamiltonian function of the 
harmonic oscillator restricted to the cotangent bundle Ao = 0, fio — reads (up to an inessential 
constant term) 

H= *h n^-A,) • (27) 

where 0^=1 (^* — A j) := f° r n = 1, (recall that Y17=i A ?/Ili=i(^ — a j) can De written 

as X^iLi^i)- Since H has the form (0), the general result proved in M allows us to infer 
immediately that 
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i) the elliptic coordinates and the corresponding momenta jij, are Darboux-Nijenhuis coor- 
dinates for a standard BH structure (Pq, Pi); 

ii) in these coordinates, P\ takes the canonical form ([!]); 

iii) the vector field X = P dH is a Pfaffian QBH vector field. The function K takes the form 

K= *h pi — w=^-w — ; (28) 

iv) a complete set of rational integrals of motion in involution is given by (|5|) . 
3.4 A solvable n— body problem 

This system belongs to a large class of integrable n-body problems in the plane, recently intro- 
duced by F. Calogero [Q. Let M = C 2n (with coordinates A := (Ai, . . . , A n ) and the conjugate 
momenta fi :— (/ii, . . . , fJ, n )) be the phase space of the dynamical system with Hamiltonian 

H ~h n^(A,-A,) ' (29) 

where gi are arbitrary smooth functions, each one depending only on the corresponding coor- 
dinate Ai, and a, b are arbitrary constants. The related Newton equations of motion take the 
form 

A fc =2V-^--6A fc . (30) 

In the case b = 0, a QBH formulation and an alternative BH formulation (according to the item 
i) of Subsec. 2.1) has been discussed in B, while in the case b = \/—lu), (uj € R) the motion 
has been proved to be completely periodic in [Q. 

Just as in the previous example of this section, by comparing (^9|) with (Q) one immediately 
concludes that: 

i) the coordinates A^ and the corresponding momenta /ii are Darboux-Nijenhuis coordinates 
for a standard BH structure (Po,Pi); 

ii) in these coordinates, Pi takes the canonical form (|l|); 

iii) the vector field X = Pq dH is a Pfaffian QBH vector field. The function K takes the form 

*-£*TEU(a^) ' (31) 
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iv) a complete set of integrals of motion in involution is given by 



Furthermore, the corresponding Hamilton-Jacobi equation is separable; a complete integral 
is S(X; bi, . . . , b n ) = —bit + W(X; b\, ... , b n ) with 




with bo = b and b\ = H . 

We wish to stress that, unlike the previous three examples, this system is not of Stackel 
type, nevertheless it is separable as well. 



4 Concluding remarks 



The Kepler and the Euler problems considered in Subsections 3.1 and 3.2 correspond to Hamil- 
tonian systems with two degrees of freedom. We recall that for such kind of systems there 
is a general result stating that a QBH formulation always exists M: however, it is essentially 
different from those presented above. Hence, the two systems are explicit examples of the non 
uniqueness of the QBH formulation. 

Indeed, in jj| one assumes to have a vector field X , Hamiltonian w.r.t. an invertible Poisson 
tensor and Liouville-integrable; passing to the action-angle variables, one can conclude that, 
for n — 2, Eq. (|J) holds for a suitable Poisson tensor P[, a function K' and an integrating 
factor p' . In particular, p' depends only on the action variables, so it is a constant of motion for 
X; this general property allows us to infer, by simple inspection, that our QBH formulations 
are different, since in both cases the functions p depend only on the configuration variables and 
therefore they are not integrals of motion for the corresponding dynamical systems. 

Finally, let us make a few comments about the relevance of the QBH formulation (we thank 
an anonymous referee for arising this question). 



Let us consider two different situations (as well as those presented in Subsect. 2.3 and Sect. 
|3[ respectively), reflecting two ways in which the QBH systems arise; these situations can be 
classified according to the fact that the phase space M of the given dynamical system X is: 

i) a BH manifold (M,Vo,Vi), with both To and V\ non invertible; 

ii) a symplectic manifold (M, P = uJq 1 ), being a symplectic tensor. 

In the first case, if one wants to solve the equations of motion through a complete integral 
of the Hamilton-Jacobi equation, then one has to pass to a symplectic manifold (as it happens 
for the stationary [^9j and the restricted flows of the KdV hierarchy, whose phase space is odd- 
dimensional). As a matter of fact, the symplectic manifold is the proper geometrical setting 
where the Hamilton-Jacobi method must be set up Q (after having been usually considered 
on cotangent bundles). To the best of our knowledge, a further generalization to Poisson (not 
symplectic) manifolds is still lacking. A possible way to achieve this goal is to perform a 
"reduction procedure" to a symplectic leaf So of one of the Poisson tensor, getting (if any) a 
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QBH formulation for X and a Nijenhuis tensor N := P\Pq (see Subsect. for details and 
notations). It is just such a tensor, living on So and not existing on M, that allows us to define 
in an intrinsic way, by means of its spectral data, the Darboux-Nijenhuis coordinates. We recall 
that such coordinates are separation variables for the Hamilton-Jacobi equation corresponding 
to any Pfaffian QBH system. 

In the second case, M is a symplectic manifold. If X is Liouville-integrable, it always 
admits (infinitely many) alternative BH formulation, not including P (as recalled in item i) of 
Subsect. 2.1). However, if one wants still to exploit Pq for getting a standard BH structure, in 
many cases a QBH formulation for X can be constructed in a quite natural way, as it has been 
shown in the examples discussed above. Of course, in this case one could also try to embed the 
given dynamical system X in a larger phase space and to lift the QBH formulation in order to 
get a BH formulation for (the lifting of) X, somehow reverting the above mentioned reduction 
procedure. However, this lifting is neither natural nor unique, and would oblige to work anew 
in a Poisson (not symplectic) manifold. 
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